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We investigate the collision cascade that is generated by a single moving particle in a static and 
homogeneous hard-sphere gas. We argue that the number of moving particles at time t grows as t^ 
and the number collisions up to time t grows as t v , with £ = 2d/(d + 2), rj = 2{d + l)/(d + 2), and 
d the spatial dimension. These growth laws are the same as those from a hydrodynamic theory for 
the shock wave emanating from an explosion. Our predictions are verified by molecular dynamics 
simulations in d = 1 and 2. For a particle incident on a static gas in a half-space, the resulting 
backsplatter ultimately contains almost all the initial energy. 

PACS numbers: 05.20.Dd: Kinetic theory, 45.50.Tn: Collisions, 47.40. Rs: Detonation waves 



The classical hard-sphere gas serves as a paradigm for 
a real molecular gas for situations where quantum ef- 
fects are negligible and where the intermolecular poten- 
tial is also negligible except when two molecules are in 
close proximity. This model has been successfully applied 
to understand many statistical properties of equilibrium 
and non-equilibrium molecular fluids 

The hard-sphere gas consists of identical spherical par- 
ticles that move along straight-line constant-velocity tra- 
jectories that are interrupted by elastic collisions with 
other spheres. Despite its simplicity and considerable 
research devoted to this model, many unexpected prop- 
erties continue to be discovered, such as the slow decay 
of velocity correlations Q and the anomalous statistics 
of collision events Q • Here we investigate a phenomenon 
of a somewhat different genre, namely, the response of a 
stationary hard-sphere gas to the perturbation of a single 
moving particle. After the first collision, two particles are 
moving and these, in turn, collide with other particles, 
ultimately leading to a remarkably symmetric cascade of 
moving particles (Fig [J). In this collision cascade, we ar- 
gue that the number of moving particles at time t grows 
as tP and the total number collisions up to time t grows 
as f, with £ = 2d/(d + 2), rj = 2(d + l)/(d + 2), and d is 
the spatial dimension. 

The above setting is just the initial "break shot" in 
an infinite billiard table. In classical billiards theory, 
a wealth of beautiful phenomena has been discovered 
about periodic and chaotic motions of a single particle 
that collides elastically with a confining boundary 0, H| . 
In real billiards, however, multiple particles move simul- 
taneously. While some intriguing results exist about 
the collision dynamics of few-particle systems @, little 
is known quantitatively about infinite-particle systems, 
such as the break shot that is the focus of this work. 

The sudden injection of energy by the incident par- 
ticle in the interior of a zero-temperature gas can also 
be regarded as an explosion, and the collision cascade is 
then analogous to the shock wave that propagates out- 
ward from the initial detonation. From hydrodynamic 
theory @, basic properties of the shock wave, such 
as its propagation velocity, and the radial dependences 
of the density and temperature profiles behind the wave 



are well understood. Strikingly, the collision cascade in 
the hard-sphere problem shares many quantitative prop- 
erties of this shock wave. 




FIG. 1: (color online) Collision cascade in two dimensions 
at volume fraction p = 0.4 and t ~ 552. Moving spheres 
are shown actual size (red) and stationary ones at half radius 
(blue) . The cascade slowly moves rightward due to a particle 
initially at the center moving to the right. 



Two natural characteristics of the cascade are N t , the 
number of moving particles at time t, and Ct, the to- 
tal number of collisions up to time t. These quantities 
depend on the initial sphere locations, but numerical sim- 
ulations indicate that relative sample-to-sample fluctua- 
tions vanish in the long time limit. Therefore we focus on 
the number of moving particles and number of collisions 
averaged over many initial conditions, N(t) = (N t ) and 
C(t) = (Ct), respectively. 

The number of moving particles can depend on time 
t, the density n, the spatial dimension d, the initial ve- 
locity vq of the moving particle, and the particle radii 
a: N(t) = N[t, Vq, a, n, d]. Since the particle number 
is dimensionless, dimensional analysis tells us that N(t) 
should depend on the three independent dimensionless 
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parameters T ■■ 



i x l d vot, p = nflda d and d: 
N = N[T,p,d], 



(1) 



where fid = n d / 2 /T(l + d/2) is the volume of unit sphere 
in d dimensions. The self-similar nature of the collision 
cascade then suggests that N(t) should have the power- 
law time dependence 



N ~ T« 



(2) 



with exponent £ that potentially depends on p and d 

We determine £ by writing a rate equation for the evo- 
lution of the number of moving particles. From Fig. [TJ 
only particles within a shell of thickness one mean-free 
path I ~ vt ~ l/(na d ~ 1 ) from the edge of the cascade at 
time t can recruit additional moving particles in a colli- 
sion time t. We now assume equipartition in which the 
initial energy is shared approximately equally among the 
moving particles; simulations clearly support this asser- 
tion. At time i, each particle in the cascade then has 
energy that is roughly TV -1 times the initial energy, and 
root-mean-square velocity v(t) ~ A -1 / 2 . Using the ob- 
servation from Fig. [1] that the collision cascade is com- 
pact, the number of particles within the collision shell is 
n£R d ~ 1 . Consequently, the rate equation is 



dN 
~~dt 



nlR a 



n J 



(d-l)/d 



(3) 



Now using v ~ N x ^ 2 and solving the resulting differen- 
tial equation yields ©, with 



2d 
d + 2 



(4) 



a universal exponent that does not depend on p. Notice 
also that N itself does not depend on p. We might an- 
ticipate that p determines the time at which asymptotic 
behavior sets in, with this crossover time becoming long 
for small p. 

As mentioned at the outset, the properties of 
the collision cascade can also be determined by a 
hydrodynamically-based dimensional analysis. The per- 
turbation caused by a single moving particle can be 
viewed as the instantaneous release of energy E = ^mv 2 . 
This point energy source is akin to detonating an explo- 
sion, and the latter problem has an elegant scaling solu- 
tion 0, [|| as long as the pressure behind the shock wave 
greatly exceeds that in front. This condition is always 
satisfied in our case, as the exterior gas has zero pres- 
sure. In this infinite Mach number limit, the radius R(t) 
of the shock wave can only depend 0] on the energy re- 
lease E, the mass density nm, where m is the particle 
mass, and the time t. Therefore the radius must scale as 



R(t) ~ (M 2 /nm) 1/(d+2) , 



(5) 



as this is the only variable combination with units of 
length. The only feature not determined by this dimen- 
sional analysis is the prefactor in Eq. (JSJ) which requires 



solving the hydrodynamic equations of motion 0, Q . Us- 
ing the result for R(t) in conjunction with E = ^mv 2 and 
N = nildR d , we again recover fl2) & (j4}. 

To obtain the number of collisions C(t) up to time t, we 
note that in a collision between a moving and a station- 
ary particle one additional particle moves Q, while in a 
collision between moving particles their number remains 
the same. Thus we have the trivial bound C(t) > N(t), 
which is valid for all times if we count the impulse that 
causes the initial particle motion as a collision. We again 
anticipate that the total number of collisions scales alge- 
braically with time C oc T" 1 . 

To determine 77, we write the rate equation for C 



dC_ 
~dt 



N 

T 



(6) 



that merely states that A additional collisions occur per 
collision time. We now use the fact that r ~ t/v is 
proportional to N 1/2 /p. Together with Eq. © for A 
and T — n^^Vot, we integrate the rate equation to find 



(7) 



Therefore 77= 1 + £/2 = 2(d+ l)/(d + 2). 

Although the initial trigger of the cascade is direc- 
tional, the bulk motion of the cascade is negligible for 
d > 2. Since the total momentum in the cascade equals 
the initial momentum mvo = 1, when the cascade encom- 
passes N particles, its center-of-mass velocity is given by 



1, or v c 



N~ 



T Thus the position of 



the center-of-mass of the cascade is: 



tf) dt' 




(8) 



To test our predictions, we perform molecular dynam- 
ics simulations on a gas of perfectly elastic hard spheres. 
In two dimensions, the system is initialized by placing, 
one-by-one, a given number of spheres of the same mass 
uniformly at random within a square box until a vol- 
ume fraction p is reached. If overlap is created when a 
sphere is introduced, one of the overlapping spheres is 
repositioned to eliminate this overlap. Without loss of 
generality we set the radius and the mass of spheres, and 
also the initial speed of the initially-moving sphere that 
is in the middle of the box to 1. The remaining spheres 
are initially at rest. We stopped the simulations before 
any moving sphere reaches the edge of the box. Our data 
for the time dependence of the average number of moving 
spheres and the number of collisions is in good agreement 
with theory in the long-time limit (Fig. [2]) . We also find 
that the exponents £ and 77 are close to the theoretical 
values of 1 and 3/2 respectively, and are independent of 
the volume fraction p for 0.1 < p < 0.45. It is worth 
emphasizing that moving particles have their directions 
quickly randomized for d > 2 so that the cascade be- 
comes symmetrical about its center of mass in spite of 
the directional initial perturbation. 
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FIG. 2: Number of moving spheres N(T) and number of colli- 
sions C(T) versus scaled time T = n 1 ^ 2 vot in two dimensions 
for volume fraction 0.2 The dashed lines have slopes 1 and 
3/2 as predicted by theory. 

The analogous problem in d = 1 is pathological be- 
cause there is no "mixing" — two elastically colliding par- 
ticles merely switch their velocities. By this construction, 
the collision cascade reduces to the initial particle freely 
traversing the particle array. To avoid this pathology, 
we generalize slightly and allow for distributed particle 
masses. This generalization leads to collisional veloc- 
ity mixing and thus represents a proper one-dimensional 
counterpart of our original system. We specifically con- 
sider uniformly distributed masses between 1 and 1 — w, 
with w > 0.2 to avoid the slow approach to asymp- 
totic behavior if the mass distribution is narrow. We 
initially place particles at integer points in the domain 
x G [— L,L]; the particle at the origin is assigned a unit 
mass and velocity vq = 1. 

With this initial condition, only particles in the re- 
gion x > are excited at early times. However, the 
cascade gradually becomes more symmetrical! After an 
initial transient, the ratio of the positions of the right- 
most and leftmost moving particles |x+/a;_| systemati- 
cally decreases and approaches a value that is close to 1 
as t — > oo. This same approach to symmetry occurs in 
the density profile of the gas. The collision cascade also 
satisfies the Rankine-Hugoniot conditions @ that relate 
discontinuities in thermodynamic properties of the gas 
across the shock front in the infinite Mach number limit. 
For example, the gas just behind the leading edge of the 
cascade is twice as dense as that of the unperturbed gas 
(for a monoatomic gas in one dimension; the density ra- 
tio is d+1 in d dimensions), in agreement with Rankine- 
Hugoniot. Similar results have been found in large-scale 
numerical simulations of a planar shock in a hard-sphere 
fluid [Ty] . Finally, as shown in Fig. [31 our simulation re- 
sults for N(t) and C(t) again agree with the theoretical 
predictions of £ = 2/3 and 77 = 4/3 in d=l. 

We now study the "splatter" when a particle is nor- 
mally incident on a stationary gas in the half space x > 
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FIG. 3: (color online) Dependence of N(T) and C(T) versus 
scaled time T = nvot in one dimension on a double logarith- 
mic scale for different widths of the distribution of particle 
masses. The asymptotic behavior is approached more slowly 
for narrower mass distributions. The dashed lines have slopes 
4/3 and 2/3 as predicted by theory. 
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FIG. 4: (color online) Collision cascade in a two-dimensional 
half space at volume fraction p = 0.4 and t — 400 when a 
particle impacts from the left at t = 0. 



(Fig. [4} . Again a collision cascade penetrates the gas in a 
manner similar to that of the bulk system, but some par- 
ticles are ultimately ejected backward. Even though this 
backsplatter is visually a small fraction of the particles 
in the cascade, the energy contained in this backsplatter 
approaches value that is close to (and it perhaps equal 
to) the initial energy. In one dimension the situation is 
more clear but still surprising. At short times the inci- 
dent particle continues moving rightward as it recollides 
with the particle array. Eventually, however, the incident 
particle is reflected so that it never experiences another 
collision. Particles 2,3,... undergo subsequent and simi- 
lar reflections so that a finite fraction of the particles are 
ejected backward as a non-interacting "fan" . The kinetic 



4 



energy in the region x > slowly decreases with time 
and seems to vanish as i _/3 with (3 « 0.1, an observation 
that can be justified heuristically Thus all the en- 
ergy ultimately resides in the region x < and, in fact, 
is concentrated in the trailing edge of the backsplatter. 
(A related issue of optimal energy transmission through 
a finite array was considered in fl2l|.) An important 
consequence of the vanishing of the energy in the region 
x > 0, is that Eq. §5§ now predicts that the radius of the 
collision cascade grows as £( 2 ~^U 3 instead of t 2 / 3 that 
arises from strict energy equipartition. 

To summarize, the initial break shot on an infinite 
billiard table with a finite density of stationary spheres 
creates an almost symmetrical collision cascade within 
which excited particles share energy nearly equally. Di- 
mensional arguments determine the extent of the cascade 
as a function of time. The properties of the cascade are 



closely analogous to those for a continuum shock wave 
that emanates from an explosion in a gas. It will be 
worthwhile to understand the limits of applicability of 
the continuum hydrodynamic analogy, especially for one 
dimension. The collisional splash problem reveals un- 
expected features, the most prominent being that ap- 
parently all of the energy is asymptotically carried by 
the backsplatter. A major challenge would be to prove 
that the energy transmitted to the initially occupied half- 
space does indeed decay to zero; if so, to then deter- 
mine the decay law and whether the splatter is described 
within a hydrodynamic framework. 
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